Entanglement distillation by dissipation and continuous quantum repeaters 
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Even though entanglement is very vulnerable to interactions with the environment, it can be 
created by purely dissipative processes. Yet, the attainable degree of entanglement is profoundly 
limited in the presence of noise sources. We show that distillation can also be realized dissipatively, 
such that a highly entanglement steady state is obtained. The schemes put forward here display 
counterintuitive phenomena, such as improved performance if noise is added to the system. We also 
show how dissipative distillation can be employed in a continuous quantum repeater architecture, 
in which the resources scale polynomially with the distance. 
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Entanglement plays a central role in applications of 
quantum information science such as quantum computa- 
tion, simulation, metrology, and communication. How- 
ever, any quantum technology is challenged by dissipa- 
tion. The interaction of the system with its environ- 
ment is regarded as a major obstacle, and in particular 
the degradation of entangled states due to dissipation is 
typically considered to be a key problem. Contrary to 
this belief, new approaches aim at utilizing dissipation 
for quantum information processes |Tj including quantum 
state engineering [2H1] , quantum computing [3] , quantum 
memories [5] , the creation of entangled states [5] , and er- 
ror correction [7]. 

Entanglement generated by dissipation has been 
demonstrated experimentally J3J following a recent theo- 
retical proposal [6] . The main advantage of this scheme 
lies in the fact that entangled states are generated in a 
steady state. Furthermore, as opposed to standard meth- 
ods, the desired state is reached independently of the ini- 
tial one. By coupling two quantum systems to a common 
environment (e.g. the the electromagnetic field S ) a ro- 
bust entangled steady state can be quickly generated and 
maintained for an arbitrary long time without the need 
for error correction such that entanglement is available 
any time. 

As any other scheme, dissipative protocols are exposed 
to noise sources, which degrade the quality of the pro- 
duced state and render it inapplicable for many impor- 
tant applications in quantum information, like quantum 
communication where noise effects increase dramatically 
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FIG. 1: (Color online) Entanglement distillation by dissipa- 
tion a) Distillation setup without communication, b) Distil- 
lation setup including classical communication. 



with the distance. By means of distillation [9], entan- 
glement can be improved at the expense of using several 
copies. In combination with teleportation, this method 
allows for the construction of quantum repeaters |10) . 
which enable the distribution of high-quality entangle- 
ment for long distance quantum communication with a 
favorable scaling of resources. Unfortunately, existing 
schemes for distillation and teleportation are incompat- 
ible with protocols generating entanglement in a steady 
state, since they require the decoupling of the system 
from the environment, such that the advantages are lost. 
Hence, new procedures which are suitable to accommo- 
date dissipative methods such that all advantages can 
be retained and used for quantum repeaters are highly 
desirable. 

We introduce and analyze different dissipatively driven 
distillation protocols, which allow for the production of 
highly entangled steady states independent of the initial 
one and present a novel quantum repeater scheme featur- 
ing the same properties. More specifically, this protocol 
continuously produces high-quality long-range entangle- 
ment. The required resources scale only polynomially 
in the distance. Once the system is operating in steady 
state, the resulting entangled link can be used for appli- 
cations. Remarkably, the time required to drive a new 
pair into a highly entangled steady state is independent 
of the length of the link such that this setup provides 
a continuous supply of long distance entanglement [10] . 
Apart from that, the proposed distillation protocols ex- 
hibit several intriguing features. We present, for exam- 
ple, a method which allows for distillation in steady state 
where none of the individual source pairs is entangled, 
and describe another one whose performance can be im- 
proved by deliberately adding noise to the system. 

In the following, we introduce two types of dissipative 
distillation protocols suitable for different situations. We 
start out by explaining scheme I which is physically mo- 
tivated and consider the situation shown in Fig. [T] Two 
parties, Alice and Bob, share two source qubit pairs Si 
and S2, which are each dissipatively driven into an en- 
tangled steady state and used as resource for creating 
a single highly entangled pair in target system T. As- 
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suming Markov dynamics, the time evolution of the den- 
sity matrix p can be described by a master equation of 
Lindblad form p = 7 [ApA^ - \ (pA^A + A^Ap)) with 
rate 7 and will be abbreviated by the short hand no- 
tation p — r )Z A {p) in the following. The entangling 
dissipative process acting on the source qubits consid- 
ered here, is the single-particle version of the collective 
dynamics realized in [8] and corresponds to the mas- 
ter equation p = £° at (p) = 7 (£ A (p) + £ s (p)) with 
A = cosh(r)<7~ licc + sinh(r)a+ ob and B = cosh(r)cr Bob + 
sinh(r)cr+ ico , where a~ = |Q> <1 j and a+ = |1)(0|. The 
unique steady state of this evolution is the pure entangled 
state \ip) = (| 00) - A|ll)) + A 2 , where A = tanh(r). 
It is subject to local cooling, heating and dephasing 
noise described by £ noise (p) = e c (& a Aiio.(p)+^Bo b (pjj + 

e h (£<.-( P )+£°'B + ob(p)) + £d (&w\u*.(p) + fii^bob^)). 

We assume that the entangling dynamics acting on Si 
and S2 is noisy, while the target system is protected (this 
assumption will be lifted below). The source qubits are 
locally coupled to T such that 

p = £«» + £°>) + £^» + + 2 Mi Jp) + £ B » , 

where £ A ii C e(£Bob) acts only on Alice's (Bob's) side. We 
choose £ A iicc(p) = — £Bob(p) = iSr\F,p], corresponding 
to the unitary evolution with respect to the Hamilto- 
nian F = ^2ij\jti s )(kjs\, where |6 B ) = |0 Sl l S2 ) and 
|l s ) = |l Sl S2 ). Note that this distillation protocol does 
not require any classical communication or pre-defined 
correlations. As can be seen in Fig. [2^l, the efficiency is 
mainly determined by the mixedness of the source states 
rather than their entanglement. In the absence of errors, 
the target system reaches a maximally entangled state. 

In order to allow also for noise acting on T, we in- 
clude now classical communication. As shown in the ap- 
pendix, any Lindblad operator of the form £ T "xc (p) = 
{Tuxc(p) — p) , where is an arbitrary LOCC channel 
[IT] , can be realized using local dissipative processes in 
combination with classical communication |12j . In par- 
ticular, this allows for the stabilization of the distillation 
schemes discussed below against errors acting on the tar- 
get system by running them using m blocks of source 
pairs, which are all coupled to the same target state (see 
Sec. 3 in [13J as illustrated in Fig. [3^. If sufficiently 
many source-blocks, to, are used, the dynamics is dom- 
inated by the desired processes. For clarity, we discuss 
the following distillation schemes in the absence of target 
errors, which corresponds exactly to the limit m — > 00. 

Thus, we consider the master equation 
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The LOOC map Tp(p) is defined by the four Kraus oper- 
ators ¥ A ®F B , ® P B , Pa ® Pb > p a ® p b = where p i p± 
are the projection onto the one excitation subspace and 
its orthogonal complement. Alice and Bob measure the 
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FIG. 2: (Color online) Dissipative distillation according to 
scheme I without communication (panel a) and including clas- 
sical communication (panels b-d). The full red lines show the 
steady state entanglement of formation (Eof) of system T ■ 
The dashed blue lines depict the steady state Eof of the source 
state si if no distillation is performed (a,c,d) and during the 
protocol (b). For better visibility the blue dashed curve is 
multiplied by a factor 30 in panels b and c. The dotted green 
lines show the entropy of si which is a measure of its mixed- 
ness. a) EoF attainable without communication versus error 
rate £n = £h = £c = £<i- b) EoF versus the noise parameter 
e c . c) EoF versus error rate £n- The black dotted curve rep- 
resents the entanglement of the total source system measured 
in log negativity, d) EoF versus the parameter r 



number of excitations on their side. After successful pro- 
jection onto the subspace with one excitation Pa <8> Pb, 
a flip operation F is performed, in the unsuccessful case 
no operation is carried out. 

As shown in Fig. [2}d, the scheme is robust against lo- 
cal noise of cooling- type (£ CT (p)). This kind of noise can 
even be used to enhance the performance of the distilla- 
tion protocol in the steady state at the cost of a lower 
convergence rate. Thus, counterintiutevly, it can be ben- 
eficial to add noise to the system in order to increase the 
distilled entanglement. Moreover, the steady state en- 
tanglement of the source pairs is zero in the absence of 
cooling noise for the parameters considered in Fig. [2]d, 
if no distillation scheme is performed. For increasing e c , 
the entanglement in Si and S2 increases, reaches an op- 
timal point an decreases again. Yet, the entanglement 
that can be distilled from these pairs is monotonously 
increasing and displays a boost effect. Panel c also hints 
at another counterintuitive effect, namely that entangle- 
ment can be distilled even though none of the source pairs 
is (individually) entangled in the steady state. This can 
be explained by noticing that the two-copy entanglement 
can be maintained for high noise rates when the single- 
copy entanglement is already vanishing. Fig. [2]i shows 
that the distilled entanglement increases considerably for 
small values of the parameter r despite the decrease in 
the entanglement of the source pairs. This is due to the 
fact that the protocol is most efficient for source states 
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FIG. 3: (Color online) Building blocks of a dissipative quan- 
tum repeater, a) Noise resistent distillation setup. The pro- 
cess acting on the target system is boosted using several copies 
of the source system, b) Continuous entanglement swapping 
procedure. 



close to pure states, where it allows one to distill quickly 
highly entangled state. 

In settings where the source states can be highly mixed, 
another distillation scheme (scheme II hereafter) is a 
method of choice and will be explained in the following. 
We analyze a generic model, which can be solved exactly 
and allows one to reduce the discussion to the essential 
features of dissipative entanglement distillation. As in 
standard distillation schemes, we study the general prob- 
lem in terms of Werner states [14] , since many situations 
can be described this way and a wide range of processes 
can be cast in this form by twirling [T^j. Werner states 
are of the simple form pw{f) — f& + (1 — /)(! — ^)/3, 
where is a projector onto the maximally entangled 
state |00) + |11), and I the identity operator. We as- 
sume a process which drives each source pair into the 
state il, p = 7 (tr(p)fi — p) = jE(p). Local depolar- 
izing noise is added in the form of the Lindblad term 
N(p) = (pAHce®S-p) + (II®pBob-p) , where p Alicc (p Bob ) 
denotes the reduced density matrix of Alice's (Bob's) 
system and 1 the normalized identity. This term de- 
scribes the continuous replacement of the initial state by 
the completely mixed one The source system reaches the 
steady state p s oc 7^ + dl of the total master equation 
p = r yE(p) + ^N(p) at least exponentially fast in 7 (see 
Sec. 4 in [13]). A continuous distillation process based on 
a standard protocol [15] can be constructed considering 
n source pairs which are independently driven into the 
steady state p s and a target system T ■ T is coupled to 
the source pairs by a dissipative dynamics of the form 
p = <5d (tr(p)To( ( o) — p), where the completely positive 
map To(p) corresponds to a process which acts on the n 
source pairs and distills a single potentially higher entan- 
gled copy. The output state is written on T, while the 
n source pairs are re-initialized in the state 1L The total 



master equation is given by 

n 

p = J2 (lEi(p) + |jVi(p)) + S D (T D (p) - p), 
i=i 

where Ei, Ni denote entangling and noise processes on 
the ith source qubit pair. The steady state has a fi- 
delity of / = f} dx foils — (fs — 0.25)x^d"), where f s and 
f D (f) = tr(nT D {p w (f)® n )) are the fidelity of p s and the 
output of the distillation protocol with n input states of 
fidelity /. High fidelities require low values of <5d. How- 
ever, the solution p{t) (see [13], Sec. 4) shows that fast 
convergence requires high values of this parameter. A low 
convergence speed on the target system is extremely dis- 
advantageous if noise is acting on T. Therefore, a boost 
of the process as illustrated in Fig. [3^, is required. This 
way, the new convergence rate is given by mSn while the 
back action on each source system remains unchanged 
(see [13], Sec. 3). 

The distribution of entanglement over large distances is 
one of the big challenges in quantum information science. 
In quantum repeater schemes, entanglement is generated 
over short distances with high accuracy and neighbor- 
ing links are connected by entanglement swapping. This 
procedure allows one to double the length of the links, 
but comes at the cost of a decrease in entanglement for 
non-maximally entangled states. Therefore a distillation 
scheme has to be applied before proceeding to the next 
stage, which consists again of entanglement swapping and 
subsequent distillation. The basic setup for a continuous 
entanglement swapping procedure is sketched in Fig. [3]d. 
It consists of three nodes operated by Alice, Bob and 
Charlie, where Alice and Bob as well as Bob and Char- 
lie share an entangled steady state. By performing a 
teleportation procedure, an entangled link is established 
between Alice and Charlie and written onto the target 
system, while the source systems are re-initialized in the 
state 1L This corresponds to LOCC operation T sw (p). 
The whole dynamics is described by the master equation 

2 

p = E (t^OO + f ^G°)) + <W(T SW ( P ) - p). 

i=l 

The steady state has a target fidelity of / = 

(27+^!x7+g 5W ) (/s™(/s) ~ 3) + b where fsw(L) is the 
output fidelity of the entanglement swapping protocol for 
two input states with fidelity f s (see [13], Sec. 5). 

The basic idea of a nested steady state quantum re- 
peater is illustrated in Figj4] At the lowest level, entan- 
gled steady states are generated over a distance L . At 
each new level, two neighboring states are connected via 
a continuous entanglement swapping procedure and sub- 
sequently written onto a target pair separated by twice 
the distance. The distillation and boost processes, that 
are required in each level to keep the fidelity constant are 
not shown in this picture. The resources required for this 
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FIG. 4: (Color online) Steady state quantum repeater scheme. 

repeater scheme can be estimated as follows. Entangle- 
ment swapping processes acting on source pairs of length 
I with fidelity // result in entangled target pairs of length 
21, with degraded fidelity fax < This reduction is due 
to the swapping procedure, noise acting on the target sys- 
tem and the back-action from entanglement distillation. 
Stabilization against noise acting on the target systems 
is achieved by coupling each of them to m copies of the 
source system and requires therefore 2m source pairs of 
length I. In order to obtain a fidelity f2i > fu n copies 
of these error stabilized links are used as input for a n 
to 1 distillation process. The distilled state is mapped 
to another target pair of length 21, which also needs to 
be stabilized against errors using m copies of the blocks 
described. Hence, in total 2m 2 n pairs of length I are re- 
quired for a repeater stage which doubles the distance 
over which entanglement is distributed. For creating a 
link of length L — L 2 k , (2m 2 n) k source pairs are needed, 
where k is the number of required iterations of the re- 
peater protocol. Therefore, the required resources scale 
polynomial with (L/ L ) Lo ^ 2m2n \ In Sec. 5 in |13j, we 
discuss a specific example scaling with (L/L ) 16A . The 
convergence time of the total system scales only loga- 
rithmically with the distance L. Once the steady state is 
reached, the entanglement of the last target system can 
be used e.g. for quantum communication or cryptogra- 
phy. The underlying source systems are not effected by 
this process and remain in the steady state. Therefore, 
the target state is restored in constant time. 

In conclusion, we have shown how entanglement can 
be distilled in a steady state and distributed over long 
distances by means of a dissipative quantum repeater 
scheme serving as stepping stone for future work aiming 
at the optimization in view of efficiency and experimental 
implementations. 

APPENDIX The continuous exchange of classical 
communication is added in the framework of dissipative 
quantum information processing, by assuming that Al- 
ice and Bob have access to a system, which is used for 
communication only and considering the master equation 




States referring to the communication system at Alice's 
and Bob's side are labelled by subscripts ca and cb- 
Alice's communication system is continuously measured 



in the computational basis yielding the quantum state 
|i CA ) with probability (i CA |p A iico|*c A ) an d reset to the state 
|0 CA ), while the communication system on Bob's side is 
set to the measurement outcome. This way, classical in- 
formation can be sent at a rate T, but no entanglement 
can be created (see [13], Sec. 2). As proven in Sec. 2 in 
|13) , any operation that can be realized by means of local 
operations and classical communication (LOCC) can be 
constructed in a continuous fashion using communication 
processes £a->b an d £b->Aj if the rate T is fast compared 
to all other relevant processes including the retardation 
due to back and forth communication. 
We thank Eugene Polzik for helpful discussions and ac- 
knowledge support from the Elite Network of Bavaria 
(ENB) project QCCC, the DFG-Forschungsgruppe 635 
and the EU projects COMPAS and QUEVADIS. 



[1] M. B. Plenio and S. F. Huelga, Phys. Rev. Lett. 88, 
197901(2002); B. Kraus and J. I. Cirac, Phys. Rev. Lett. 
92, 013602 (2004); F. Benatti, R. Floreanini, and U. Mar- 
zolino, Phys. Rev. A 81,012105 (2010); R. Bloomer, M. 
Pysher and O. Pfister, |arXiv:1007 .2369 (2010); J. Cho, 
S. Bose and M. S. Kim, Phys. Rev. Lett. 106, 020504 
(2011) (2010); M. J. Kastoryano, F. Reiter, and A. S. 
S0rensen, Phys. Rev. Lett. 106, 090502 (2011); A. Mari, 
J. Eisert, |arXiv:1104.0260"] (2011); J. T. Barreiro et al., 
Nature 470, 486 (2011). 

[2] J. F. Poyatos, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 
77, 4728 (1996). 

[3] S. Diehl et al. Nature Phys. 4, 878 (2008). 

[4] F. Verstraete, M. M. Wolf, and J. I. Cirac, Nature Phys. 
5, 633 (2009). 

[5] F. Pastawski, L. Clemente, and J. I. Cirac, Phys. Rev. A 

83, 012304 (2011). 
[6] C. A. Muschik, E. S. Polzik, and J. I. Cirac, 

|arXiv: 1007.22091 (2010), see also: A. S. Parkins, E. 

Solano, and J. I. Cirac, Phys. Rev. Lett. 96, 053602 

(2006). 

[7] Joseph KerckhofF, Hendra I. Nurdin, Dmitri S. Pavlichin, 
Hideo Mabuchi, Phys. Rev. Lett. 105, 040502 (2010); J. 
P. Paz and W. H. Zurek, Proc. of the Royal Soc. of Lon- 
don. Series A 454, 355 (Jan. 1998). 
[8] H. Krauter et al. |arXiv:1006.4344| (2010). 
[9] C. H. Bennett et al. Phys. Rev. Lett. 76 722 (1996). 

[10] H.-J. Briegel, W. Diir, J. I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 81, 5932 (1998). 

[11] LOCC channels are completely positive trace preserving 
maps that can be realized by means of Local Operations 
and Classical Communication. 

[12] We assume, that the time scales for classical communi- 
cation T" 1 (see appendix) are sufficiently long such that 
retardation effects can be ignored. 

[13] Supplemental Material 

[14] R. F. Werner, Phys. Rev. A 40, 4277 (1989); M. 
Horodecki and P. Horodecki, Rhys. Rev. A 59, 4206 
(1999). 

[15] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. 
K. Wootters, Phys. Rev. A 54, 3824 (1996). 



5 



SUPPLEMENTAL MATERIAL 

Here, we explain the results presented in the main text. 
In Sec. 1 and Sec. 2, two dissipative distillation schemes 
are discussed. Scheme I is suited for settings where a dis- 
sipative processes is available which produces entangled 
steady states, that are close to pure states. If only very 
mixed steady states are available as input, scheme II is 
preferable (which we explain in detail in Sec. 4). In Sec. 2, 
the notion of continuous exchange of classical information 
between two parties is introduced in the master equation 
formalism and it is shown that arbitrary LOCC channels 
can be realized using local dissipation and classical com- 
munication. In Sec. 3, we explain how the continuous 
protocols used here can be made robust against noise. 
Finally, in Sec. 5, we analyze the dissipative quantum 
repeater scheme put forward in the main text in detail. 



1. SCHEME I: DISSIPATIVE ENTANGLEMENT 
DISTILLATION FOR SOURCE STATES CLOSE 
TO PURE STATES 

In this section, we explain two variants of scheme I (STj. 
In Sec. 1.1, we discuss a protocol, which allows for dissi- 
pative entanglement distillation without communication. 
Sec. 1.2 is concerned with a related protocol, which in- 
cludes classical communication. Both protocols produce 
Bell-diagonal steady states which can further distilled us- 
ing scheme II presented in Sec. 4. 



1.1 Dissipative entanglement distillation without 
communicat ion 



We consider the setup illustrated in Fig. |S.1| The dis- 
sipative dynamics driving the two systems si and S2 is 
physically motivated and can be implemented by cou- 
pling the systems located at Alice's and Bob's side to a 
common bath, for example the vacuum modes of the elec- 
tromagnetic field J5J |S] . The entanglement which can be 
attained per single copy is limited for a given dissipative 
process. Moreover these systems are subject to noise. 
Still, it is possible to use these two copies as resource for 
creating a single highly entangled pair in target system 
T ■ In the absence of undesired processes, the dynamics 
described by the master equation p = 7 (£ A (p) + £ s (p)) 
(see main text) drives systems si and S2 into the state 



! cx|00y00)^A[|00)Jll) a2+ |ll) 
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Alice and Bob share a maximally entangled state |4'o) = 
(100)^11)^+111)^100)^) /V2 in a subspace with one exci- 
tation on each side. Scheme I is based on the extraction 
of entanglement from this subspace and its subsequent 
transfer to the target system by means of the flip op- 



as 



sh m 
'3 1 

Oh 

o 



o 



1 (£x.»+£ fl . 1 M) 



FIG. S.l: (Color online) Scheme I, dissipative entanglement 
distillation without communication. 



I Is) = |lsiO S2 ). Systems Si and S2 are permanently driven 
back to an entangled state. In contrast to standard distil- 
lation protocols for pure states |S10j , the presence of this 
strong process leads to a substantial decrease in the en- 
tanglement if the flip operations on Alice's and Bob's side 
are not applied simultaneously. Hence, the coordination 
of their actions, e.g., using fast classical communication, 
seems to be essential. Surprisingly, the desired dynam- 
ics can be realized in the absence of communication or 
predefined correlations using local unitary evolutions. 

This is possible by exploiting the symmetry of the max- 
imally entangled state l^o). More specifically, \^q) is in- 
variant under any unitary operation of the form U ®U , 
while less entangled pure states are not. U denotes the 
complex conjugate of U . Such an operation can be im- 
plemented without communication as the time evolution 
of a sum of local Hamiltonians H = Ha <8 I — I ® Hb- 
Here, we use the flip operation such that the correspond- 
ing master equation is given by 

p = 7 (£ A ^ (p)+£ B - (p) + £^ (p)+£ B - (p)) 
+ iS v [F ® I - 1 ® F, p] 
+ e c (£ a ^ (p) + £ fcs i (p)+£ a ^ (,o)+£^ (p)) 

+ s h (£< (p) +£< (p) +£< (p) (p)) 
+ e d (£<^(p) + £ 6 4^(p)) , 



eration F = |jr*s)(*rjs|, where 



S = 



o sl i S2 ; 



and 



where a — (T A ^. co and b = cr B ^ jb . The first line corresponds 
to the entangling dissipative process (described by 
nonlocal jump operators A and B) acting on the two 
source systems as explained in the main text. The 
second line describes the unitary coupling of the target 
system to the entangled subspace of the two source 
systems and the last three lines represent undesired 
processes. More specifically, we include dephasing at a 
rate ed as well as noise terms, which create (annihilate) 
excitations locally at the heating (cooling) rate eh (e c ). 
Note that the noise types considered here also include 
depolarizing noise. The target system itself is assumed 
to be protected (below, a variant of this scheme is 
described, which includes classical communication and 
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FIG. S.2: (Color online) Scheme I, dissipative entanglement 
distillation including classical communication channels. 



can be made robust against noise acting on the target 
system). 

A disadvantage of the unitary evolution employed 
here lies in the fact that the source system is subject 
to a back-action of the target state, which depends on 
the quantum state of T. Accordingly, the evolution of 
the source systems is highly dependent on the state of 
the target pair. It remains an open question, whether 
schemes, similar to the one described in Sec. 3 can be 
used to render this protocol robust against errors on 
the target system, or whether this is a special feature of 
protocols including classical communication. 



1.2 Distillation using scheme I including classical 
communicat ion 



We consider the setup illustrated in Fig. |S.2| As ex- 
plained in Sec. 1.1, the dissipative entangling process act- 
ing on the source systems s± and S2 has the property that 
Alice and Bob share a maximally entangled state if the 
resulting steady state is projected onto the subspace with 
one excitation on each side. Ideally, this quantum state is 
then transferred to the target system T by means of the 
flip operation defined above. In this subsection, we in- 
troduce a classical communication channel, which allows 
Alice and Bob to coordinate their actions such that flip 
operations on both sides can be performed in a synchro- 
nized fashion if both sides have successfully accomplished 
a projection onto the relevant subspace with one excita- 
tion. As explained in Sec. 2, Lindblad terms of the form 
£ Tlocc (p) = {Tujcc(p) - p), where T lJXC is an arbitrary 
LOCC channel [IT], can be realized by means of local 
dissipative processes and classical communication. 
As explained in Sec. 3, this protocol is resistent against 
target errors if it is coupled to sufficiently many blocks 
source pairs. For simplicity, we explain here the basic 
protocol in the absence of target errors, which corre- 
sponds to the limit of using infinitely many source blocks 
(entanglement distillation for a finite number of source 



blocks and finite error rates is analyzed in Sec. 3 and 
Sec. 5). Classical communication allows for the imple- 
mentation of the scheme outlined above. The LOCC dis- 
tillation operation corresponding to this process, Tp, is 
given by 

T r (j>) = F A ®F B pF A ®F B 



+Pa< 



P^ P P A ® Pb 



PspPi < 

PbPPa 



Pb 

>Pb, 



where P = |0 Sl l S3 )(0 Sl l S2 | + |l Sl 0, 2 )(l Sl S2 | is the pro- 
jector onto the subspace with one excitation, and P^- = 
1 — P the projector onto the subspace with zero or two 
excitations. Note that only the first term has an effect 
on the target system. The flip operation leads to a back- 
action on the source system, which depends on the state 
of T. In order to simplify the discussion in Sec. 3, we in- 
troduce here a slightly modified version of this protocol, 
Tf(p), which does not exhibit a state-dependent back- 
action. This can be avoided by applying a twirl P3] on 
the target system prior to the flip operation 



If(p) 



1 



3 

^ 16 

ij=0 

+P A ® 



F A ® WbUhpULWa ® F 



PbPPa 

PspPi < 



>Pb 
Pb 



-Pi ® P&pP% ® P 1 



B ■ 



Uij = erf <£) a? is a unitary operation acting on the tar- 
get system only, where <ii denote the four Pauli matrices 
and <7o is the identity. Due to the twirl, this protocol 
features an enhanced back-action on the source, which is 
independent from the target. It turns out that the perfor- 
mance of this protocol is qualitatively the same as shown 
in Fig. 2 in the main text. The total master equation is 
then given by 

p = 7 (p)+£ B - ( P ) + (p)+£ B ^ (p)) 
+ 5 ¥ (T'(p) - p) 

+ e c (£ a ^(p) + £ b ^{p)+£ a ^(p)+£ b ^{p)) 
+ s h (£< (p) +£< (p) + 2< (p) + 2 bl * (p)) 

+ £d (£< a ^(p)+£ fc ^(p)), 

where a = cr.T. „ and b = <xT, . 



2. CLASSICAL DISSIPATIVE CHANNELS AND 
DISSIPATIVE LOCC 

Classical channels are easier to realize experimentally 
than their quantum counterparts and can for example be 



7 



% CB 



FIG. S.3: (Color online) Realization of a classical dissipative 
channel. 



implemented using optical fibers. Since classical chan- 
nels are insufficient for the generation of quantum cor- 
relations, long-range links can be established over large 
distances using the toolkit of classical error-correction. 
The class of LOCC operations, i.e. quantum operations 
that can be performed using local operations and classi- 
cal communication, is of essential importance in quantum 
information theory, especially in the context of entangle- 
ment distillation protocols. 

In this section, we introduce the notion of classical 
channels in the framework of dissipative quantum infor- 
mation processing. This allows us to formulate gener- 
alized LOCC operations in a continuous dissipative set- 
ting, which includes a wide range of continuous distilla- 
tion protocols. 



2.1 Classical dissipative channels 

We start out by introducing a dissipative classical com- 
munication channel. Both parties, Alice and Bob, each 
have access to a d-dimensional system which is used ex- 



clusively for classical communication (see. Fig. S.3). The 
master equation 




I PAlice I ic A ) 1 CA i CB ) (0 CA i CB | - 



=rc A -^B(p). 



describes a one-way classical communication channel. 
States referring to the communication system at Alice's 
and Bob's side are labelled by subscripts Ca and Cb re- 
spectively. Alice's communication system is continuously 
measured in the computational basis yielding the quan- 
tum state \i CA ) with probability (i CA |/0 Alicc |i CA ) and reset 
to the state |0 CA ), while the communication system on 
Bob's side is set to the measurement outcome. This pro- 
cess can be written in the form 

p = r£ A ^ B (p) = r(T(p)-p), 

where the completely positive map T(p) is an entangle- 
ment breaking operation [S2j , which maps any state to a 
separable one. The solution of this master equation p(t) 
is given by 



pit) = p(0)e 



-rt 



drT( j o(r))e r(T -* ) . 



The second term is separable, since T{p) is entanglement 
breaking. Accordingly, the classical channel introduced 
above does not produce entanglement. Moreover any en- 
tanglement present in the state p(0) is exponentially sup- 
pressed. 



2.2 Generation of Lindblad operators of the form 

T(P) - P 

In the following we prove that any dissipative time 
evolution which satisfies a master equation of the form 
p = j(T(p) — p) can be designed by means of local dis- 
sipative processes in combination with the classical com- 
munication channels introduced above in the limit of high 

Alice 



rates T. The basic setup is sketched in Fig. S.4 



and Bob hold a bipartite system, which we refer to as 
the main system. In addition both parties have access to 
several classical communication channels and can apply 
dissipative dynamics acting on the classical channels and 
their part of the main system. This setting allows for a 
wider class of dissipative evolutions on the main system 
which includes dissipative LOCC processes. In particular 
we state the following: 

Let T(p) be any LOCC map. Let £(/?) be any bounded 
Lindblad operator, i.e., max p ||£(p)|| = 1, acting on the 
main system at a rate 7. Let Alice and Bob have access to 
classical communication channels as described above. If 
both parties can apply any dissipative process of Lindblad 
form on their side, an effective dissipative time evolution 
on the main system satisfying the master equation 



p = 1 £(p) + S(T'(p)-p) 



(1) 



separable 



after an initial waiting time of the order y can be realized. 
The completely positive operator T / (p)=T(p)+0(y/a) is 
an imperfect realization of T(p) up to an error 0(y/a), 
which vanishes for small a — %r , where 7' = 7+<5 . 0(f(a)) 
denotes any hermitian (time and state dependent) opera- 
tor with a trace norm scaling with f(a) in the limit a— >0. 

Since max p ||£(p)|| = 1, the strength of the process is 
completely encoded in 7. £(p) can include a dissipative 
LOOC map itself, as discussed at the end of this sec- 
tion. The error 0(y/a) of the LOCC map is small for 
a«l. A time evolution satisfying Eq. ([!]) can be ei- 
ther obtained by starting from certain initial conditions, 
or after an initial waiting time on the order of j, dur- 
ing which no external control is required. If T ^> j + 5 
(a<ti 1), the system evolves approximately according to 
p= 7 Z(p)+6(T(p)-p). 

Note, that LOOC operations are extremely hard to pa- 
rameterize. It is known, that they can be written as a 
separable superoperator T(p) — ^ i A i ® BipA\ ® b\, 
but not every separable superoperator is a LOCC map. 
Practically, a general LOCC map can only be charac- 
terized by fixing the number of communication rounds 
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FIG. S.4: (Color online) Using local dissipative processes and 
fast classical communication, arbitrary LOCC channels can 
be implemented in a continuous fashion. 



between Alice and Bob and to specifying the exact op- 
erations that Alice and Bob perform in each round. The 
most general operation Alice and Bob can apply is a pos- 
itive operator valued measurement (POVM). This covers 
any completely positive map as well as measurements, 
unitary evolutions, etc. A POVM is specified by a num- 
ber of Kraus operators Ai, corresponding to the possible 
measurement outcomes i, where the normalization con- 
dition J^i A\Ai = I guaranties that the probabilities for 
the different possible outcomes add up to 1. 

We consider the following situation. Alice performs 
a first POVM Ai and sends her result i to Bob. Bob 
chooses a POVM depending on Alice's result i. Sub- 
sequently, he sends the result j to Alice, who chooses her 
next POVM A % j! which may depend on all previous mea- 
surement result. This procedure can be repeated many 
times. This corresponds to the application of the opera- 
tion 



E 



*D,ia,*iiJii 



,VJr 



..PX- 



t 

iOjJO,*l,3l,".)*r >3r> 



where each Kraus operator is of the form 



X, 



«0 JO:H Jli-.VJr 



__ £J«0j0,»l — «n 



and represent one possible set of measurements outcomes 
«o, JO) 2*i>ii) • • • > Vj jr for all r POVM measurements. 
Due to this lack of a concise notation for a general 
LOCC map, a complete proof of this statement would 
be lost in notation and it would be hard for the 
reader to understand the main idea of the proof. We 
restrict ourselves therefore to LOCC maps with one 
communication round, i.e., Alice sends one message 
to Bob and Bob can send an answer back to Al- 
ice once. A generalization of the following proof to 
a LOCC map with a finite number of communication 
rounds m is straight forward and will be discussed below. 

Let T(p) denote a LOCC map with one round of 
communication. This map can be realized in the 
following way: 



• Alice applies a POVM measurement with Kraus 
operators A^, obtains the measurement result i and 
sends it to Bob. 

• Bob performs a POVM measurement B^, which can 
depend on i, and sends the result j to Alice. Since 
we assume that Alice is memoryless, Bob also sends 
the measurement outcome of Alice's measurement 



• In the last step, Alice can apply any completely 
positive map on her side. This map can depend 
on both, i and j. 

Note, that Tij(p) = J2 k C kP( C kV is also a POVM map 
with Kraus operators C % u , where the measurement re- 
sults k are not used. Let Alice and Bob have n different 
measurement results for each POVM, where n can be up- 
per bounded by the square of the dimension of the sys- 
tem. For typical distillation protocols on qubits, n = 2. 
Note that all indices for Kraus operators run from 1 to n, 
and do not start with 0. This choice allows for a shorter 
notation later on (the index is reserved for indicating 
that the classical channel is operable) . The basic setup is 



illustrated in Fig. S.4 Alice and Bob have access to clas- 
sical one-way communication channels labelled CI and 
C2. CI and C2 can be used to send information from Al- 
ice to Bob and vice versa respectively. Apart from these 
classical channels, Alice and Bob hold a system subject to 
a dissipative evolution described by the Lindblad opera- 
tor £(p). In the following, this system is referred to as the 
main system. The first classical channel needs to store 
all possible measurement outcomes obtained by Alice, 
whereas the second one needs to store the measurement 
results obtained by both, Alice and Bob. We assume 
therefore that CI and C2 are n + 1 and n' + 1 = n? + 1 
dimensional systems respectively. Note, that the state 
|0) will be used to indicate that the channel input or 
output is "empty", while the states . . . , |n) represent 
n possible measurement results of Alice and the states 
|1), . . . , |n') encode the n 2 different measurement results 
obtained by Alice and Bob. The corresponding master 
equation is given by 

p = 7< c(p) + r£ A ^B(p) + r£ A ^B(p), 

where 7£(p) ^ s a process acting on the main system only. 
We assume, that the time scales for classical communi- 
cation r _1 arc sufficiently long such that retardation ef- 
fects can be ignored. The four systems used for classical 
communication are denoted by I a ,Ib,O ai Ob as shown in 
Fig. |5^1 / and O stand for "Input" and "Output". 

As a next step, local Lindblad operators are added, 
which correspond to the application of a LOCC map de- 
pending on the registers of the classical channels. The 
following three terms are added, one for each step of the 
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protocol outlined above. The first term is given by, 

S J2 £ AMm '°(p), 

z=l,fc=0 

where acts on Alice's part of the main system and 
|i)(0U on Alice's side of the first classical system, i.e., 
the input system of the first classical channel. Note, 
that we use here the short hand notation £ A (p) = 
7 (ApA^ - | (pA^A + A^Ap)) that was already intro- 
duced in the main text. This corresponds to the first 
step of the realization of the LOCC map. Alice performs 
a POVM and writes the measurement result onto the in- 
put system of the classical channel. As second term, we 
add the Lindblad operator 

ji= l,xy— 

where Bj acts on Bob's part of the main system, |0)(i| Ot 
on the output of the first classical channel and \ j, i) (x, y\ Ib 
on the input of the second channel. Note, that the sec- 
ond channel can store both values i and j at the same 
time. \j, i) stands for any encoding of i,j in the n 2 + 1 
dimensional state space, where the label zero is reserved 
for indicating the status of the channel. The summation 
over x, y starts from zero, i.e., includes the reserved zeros 
term as well as the n 2 possible measurement results. Bob 
only carries out a POVM measurement, if he receives the 
message i via CI. Afterwards he writes i, j onto the clas- 
sical channel. Note that the sum over xy implies that 
Bob overwrites any previous state of the classical com- 
munication system. The last term to be added is given 
by 

jik — l 

where acts on Alice's quantum system with T^(p) = 

J2k @k P@k ana - |0) (*) j\o ac t on ou tP u t 01 the sec- 
ond classical channel. Alice receives the message ij and 
reacts by applying T^- to complete the LOCC map. The 
sum starts from ij = 1, i.e., Alice acts only if a message 
has arrived. She does not act if the register is empty 
(|0)). Hence, the total master equation is given by 

p = 7 £(p) + rC A ->B + rC A <-B (2) 

i=l,fc=0 

+r V £ B ^l )vj| o& ®li>CKzyL( p ) 

ji=l,xy—0 

+r £ c ^®l°^ i,J 'l°-(p) 

jik—l 

The basic idea can be described as follows. The term in 
the second line starts the process of realizing T(p) at a 



rate 8 (Alice performs the first step). The following steps 
are performed with a high rate T, such that the state of 
the quantum system stays approximately constant dur- 
ing the time needed to complete the whole operation. So, 
practically, the whole LOCC map T(p) is applied at once 
at a rate S. 

In the following, this will be proven rigorously by con- 
sidering the effective evolution of the main systems after 
tracing out the classical channels. The reduced state of 
the main system can be written as 

Pm = ^ Pv kl 

ijkl 

with pij k i = (ii a jo b ki b lo a \p\iiJo b kiJo a ), where i, j 
(k, I) denote the computational bases for CI (C2). 
The indices are arranged such that their order corre- 
sponds to the order in the communication cycle, i 
refers to the input of Alice's side, j to the output 
on Bob's side, k to the input on Bob's side and / to 
the output on Alice's side. A system of differential 
equations for all piju can be derived using pijki — 
( i i a jo b ki b lo a \p\ii a jo b kiJo a ) and Eq. The desired 

terms p ooo, Pmoo, Poioo, Poo(y)o, Pooofe) evolve according 
to 

PO0O0 =7"£(POOOo) — <!>AXX)0 + T Ty (p000(x2/))j (3) 

xy—\ 

PiOoo = 7£(Piooo) - (r + 6)p i0 Q +5}^ A^k^A] 

k=0 

+ r Y T xy(PiOO(xy)), (4) 
xy=l 

Poioo = j£{poioo) - (r + <5)po;oo + r Y PtkOO 

fe=0 

+ r ^2 T *y(P0i0(xy)), (5) 
xy=l 

Poo{ij)o = 7-2(Poo(<j)o) - ( r + ?>)p00(ij)0 (6) 
+ Tj2B]po t{xy) oB^+T 

T xy(P00{ij){xy)), 

xy—0 xy—0 

PO0O(ijQ=7^(P0O0(ij))~( r ^)P0O0(i3)+ r X^ 00 fa')( x g> ^ 

xy—0 

All other terms correspond to small errors. In the first 
step of the proof it is shown that after an initial wait- 
ing time, only the states poooo, Pmoa, Pomo, Poo(ij)o, Pooo(ij) 
are significantly populated, while the population of all 
other states is small. In the next step it is shown that 
p M w poooo- In the following we will use the short-hand 
notation p := poooo- 

Bounds for occupation probabilities 

Let us define the probabilities p^ki = ^{pijkl)- A sys- 
tem of differential equations Pijki = to{pijkl) for these 
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probabilities can be derived by from the differential equa- 
tions for pijki- Traceless terms such as j£(p) and do 
no longer appear. 

We define Pooox, • ■ ■ ,Pxxxo,Pxxxx as the sum of 
Pijki, where all indices marked with X are summed 
from 1 to n, n', in order to remove the dependence 
on the Kraus operators by virtue of their normaliza- 
tion condition (X^iLi^l^i = I)- A successful ap- 
plication of the LOCC map corresponds to the se- 
ries Pxooo, Pox oo,Pooxo, Pooox, Poooo (similarly to a X- 
excitation which created on the first index, travels to 
the right and disappears in the end), poooo takes high 
values, while the other probabilities are on the order of 
which indicates that this process is fast. However, 
the situation considered here does not correspond to this 
ideal case because Alice can start a new round before 
the last one is finished, which gives rise to probabilities 
which are denoted by indices with more than one X, e.g. 
Pxxoo,Pxoxo, ■ ■ ■ which results in an incorrect realiza- 
tion of the LOCC map. Since we are not interested in the 
complete solution but only in upper and lower bounds, 
we further simplify the system by defining the probabil- 
ities Poooo,Pxbee, Poxes, Pooxh and pooox, which cover 
all possible events. X indicates that the corresponding 
index is different from zero. Therefore, the summation 
runs from 1 to n, n' . 2 stands for an arbitrary value, i.e., 
the summation starts from zero. These five quantities in- 
clude also non-ideal process, with two or more X entries 
which correspond an errors and evolve according to 

Poooo = ~ ^Poooo + Ipooox, 
Pooox = —(5 + r)pooox + TpooXH, 
Pooxs =-(S + T)p 00X E + Tpoxss, 

POXEH = —(5 + T)p QX EE + Tpxz.EE, 

Pxeee = -Fpxeee+(Kpoooo+Pooox+Pooxe+Poxee)- 
The solution shows that the steady state (ss) with 



ST 3 



5T 2 



Poooo 



POXEE 



(r+5) 4 ' 

ST 
(T+f' 



Pooox~ 



Pxe 



5 



P00XE z 



(r + <*)3' 

(8) 



is reached up to an error smaller than 0(a 2 ) after a time 
of the order of | [S3]. In the steady state, p^ooo = 1 — 
4^+0(^). Next, bounds for pxooo,Poxoo,Pooxo, Pooox 
are derived. According to Eq. Q, pxooo = -Ipxooo + 
<5poooo + rpxoox ■ Assuming that Poooo i s reached after a 
time t' , 

PxoooW = e- r( "Wooo(*0 + J t dTeriT ^hf ^ + r Pxoox 



ST 3 



where h(t,t') > is a positive function since pxoox > 
0. Hence, doubling the initial waiting time guarantees 



Pxooo 



> 



(5r J 



According to Eq. (5), poxoo = ~(T 



S)poxoo+^Pxooo+^Pxxoo+^Poxox ■ By integration, using 
the bound for pxooo and assuming that the contributions 
from poxox and pxxoo sum up to a positive function we 
conclude that rp^js < Poxoo is fulfilled after waiting for 
another period on the order of 1. Similarly, one obtains 



(r-f^e — Pooxo and 



< 



Pooox ■ Since p ooo + Pxooo + 



Poxoo+pooxo+Pooox > l- 1 ^", any probability with more 

*v 12<5 2 

r 2 



r+<5) 
>1- 

than two X entries is smaller than ^f-. In summary, 



S n b 2 6 
p — ' — < pxooo, Poxoo, Pooxo, Pooox < p, 

after a time of the order | , where the upper bounds are 
found using Eq. Hence, a steady state is reached 

where states labelled with one (more than one) X are 
occupied with probability 0(a) (0(a 2 )). 



Differential equation for po 

The evolution of p = Poooo is governed by Eq. (??). 
After a period of the order i, po = -f'(D(l + a) since 
1 1 Pooox 1 1 = Y-+0(a 2 ). Hence, for a«l, p is approxi- 
mately constant on time scales that are short compared 
to 7' . In order to obtain an equation which depends only 
on po, we solve successively the differential equations for 
ftooo, Ahoo,Poofe)o and Pooo(tf)- According to Eq. Q, 

Piooo = -rpioao+o~Aip Q A\+% 

where 91 = 7-C(p i0 oo) - <5piooo + S J^k A iPkoooA\ + 
^J2xy=i T xy(Pioo(xy)), which is bounded by j'0(a). For 
the first three terms we use pxooo — 0(a), the last term 
can be bounded by pxoox — 0(a 2 ) such that 



/0iooo(*)=/>iooo(O)e" 



-rt 



dre T{ 



r -*) (5A iPo Al 



01 



The integral /„ W^-^Ol can be bounded by 0(a 2 ). 
The initial term is suppressed by e~ Tt and therefore 
smaller than 0(a 2 ) after the initial waiting time. Hence, 



(t) = 0(a 2 ) + S / dre v ^A iPo A\. (9) 



PiOOOK 



Since the integral is mainly determined by terms close to 
r = t and po varies little on small time intervals, po can 
be assumed to be constant. To prove this, we consider 



T,t)dT, (10) 



/ X{T,t)dT= [ X(T,t)dT+ [ X( 

Jo Jo Jf 

with X(T 1 t) = e r( - T - t *>A l poA\ and t' = t- 
Since -=7 is the typical time during which po changes, it 
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is nearly constant during the interval (t,t'). From po = 
0(~f'(l + a)), we obtain that for any f" G [*', t] 

Po(0 = Po(t) + I drpoir) = p (t) + O(y/o). (11) 



The integral from to t' in Eq. arc suppressed 

at least by a factor e~VT < a (A iPo Al is on the or- 
der of one) . Inserting E q. ([IT] ) in Eq. J9l and using 
5f* dTe r< - T -^ = f (1 - e-^ 7 ^) with e - v ^ < a yields 

PioooO) = p A iPQ {t)A\ + -O(v^), (12) 



which shows that for small a, Alice applies her first 
POVM with high accuracy. 

Next, we consider the evolution of poioo ( see Eq- ([sj ) , 



POiOO 



(13) 



where 01 can be bounded by 'y'O (a) using poxoo = 
0(a), pxxoo = £>(a 2 ) and p xox = 0{a 2 ). In- 
serting Eq. Q yields p Ql0 o(t) = AM>o(0)e~ rt + 
J^Te r( - T ^(5A i p (T)A\+SO(y^)+m^ . As before, the in- 
tegral over 01 and the first term can be bounded by 0(a 2 ) 
after a waiting time. If the remaining integral is split as 



in Eq. (10), we obtain one part, where po is nearly con- 



stant and one vanishing part. The main error is again 



due to expression ( 11 ), leading to 



S S 
Po,oo(t) = -A iPQ {t)A\ + -0{yfa). (14) 



Hence, for small a, sending classical information to Bob 
causes only marginal errors on the main system. 
Next, we consider the evolution of Poo(y)oo (Eq. 

Poo(ij)o = —Fpoo(ij)a + rBjp om oBf + 01, 
where 01 can be bounded by j'0(a) such that 

Poofe)o(t) = 6 f B)A iPQ {t)A\Bf + ^O{y/o), 

which corresponds to a process, where Bob applies his 
part of the POVM and writes i, j onto his classical input 
register. Similarly, Eq. ([7]) leads to 

S S 
Ax»(«)(t) - f B t J A lPo (t)A\Bf + -O(VS), 

which corresponds to a transfer of the classical measure- 
ment results i,j back onto Alice's side. Finally, these 
results can be applied for calculating p , 

Po= n£(po)-8po+T}T ij (pooo(ij))= 7£(pp) - 5 (T' (p ) - p ) . 

ij 



m5£ s (p T ) 




FIG. S.5: (Color online) Stabilization of dissipative protocols 
against noise acting on the target system by coupling several 
source system to the same target. 



T'(p) = T(po) + 0(y/a) represents a noisy version of the 
desired LOCC map T(p). The undesired contribution 
can be suppressed by choosing a small, i.e. by choosing T 
large enough. The generalization to more than one round 
of communication is straight forward. By summing over 
the indices of the corresponding Kraus operators, one 
obtains equations for the probabilities, which are inde- 
pendent of the POVMs applied in the protocol. From 
these equations it can be concluded that only the rele- 
vant state responsible for the application of the LOCC is 
populated, while all others are suppressed by a factor of 
order a 2 after a initial waiting time. By successive inte- 
gration as shown above, the desired approximation for p 
is obtained. Many distillation protocols only require only 
a small number of rounds to reach high fidelities and often 
even only one-way communication (half a round) [9,15]. 



3. STABILIZATION OF DISSIPATIVE 
DISTILLATION SCHEMES AGAINST ERRORS 
ACTING ON THE TARGET SYSTEM 

In this section, we explain how the distillation schemes 
presented in Sec. 1.2 and Sec. 4 can be made robust 
against noise acting on the target system. The same 
method for stabilization against errors is applicable for 
both protocols and a wide range of other dissipative 
schemes, which include classical communication. The 
basic idea is illustrated in Fig. |S.5| A dissipative pro- 
tocol is run using m blocks of source systems in parallel, 
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which are all coupled individually to the same target sys- 
tem. This way, a boost effect on the desired dynamics of 
the target system can be achieved, while the back-action 
on the source pairs remains unchanged. If sufficiently 
many source systems are provided, the dynamics on the 
target system is dominated completely by the desired dy- 
namics. In the following, we explain the application of 
this method first for schemes of the type described in 
Sec. 4 and 5 and discuss then briefly the stabilization of 
scheme I. We start out by considering a target system T 
and a source block consisting of n pairs. An entangling 
dissipative process described by the Lindblad operator 
5£(p) acts on each source pair separately such that each 
of them is individually driven into an entangled steady 
state p s . The effective master equation for the target sys- 
tem, which is obtained by tracing out the source system, 
is given by 

PT = $£s(pt), 

where the Lindblad operator £s(pt) may be time 
dependent. It does not depend on the state of T but 
only on the state of the source system as indicated by 
the subscript S. Accordingly, the convergence speed 
at which £s{pt) converges to a constant operator is 
given by the rate at which the source system reaches a 
steady state. The convergence rate of the source system 
is limited by the rate 5 at which the flip operation 
mapping the quantum states of the source system to T 
is performed (see Sec. 1 and Sec. 4). 

We assume that m identical source system S\ , . . . , S m 
are individually coupled to a single source system T 
through the Lindblad operator 8 Y^iLi £(PT,Si), where 
£(pr,Si ) is a Lindblad operator acting on T and the ith 
source system (see Fig. S.5 for a schematic overview). 
We assume that these operators are identical £$. = 2 S 
such that the dynamics of the target system is governed 
by the reduced master equation 

PT = S^2 £ Si (p T ) = m8&s (pr) ■ 



each swap operation in a standard state, for example 
the identity (strict equality requires in principle also 
that all source systems start from the same initial state. 
However, different initial states have only an effect on 
the time evolution in the beginning. The following 
discussions are only concerned with the steady state of 
the system, which is independent of the initial condi- 
tions). Scheme I including classical communication (see 
Sec. 1.2) exhibits a weak state dependent back-action. 
As explained in the end of Sec. 1.2, this can be avoided 
by applying a twirl [T3] on the target system prior to 
each flip operation. Hence, the stabilization method 
outlined above is directly applicable to this modified 
version of the scheme |S4j . 

By boosting the desired dynamics on the target 
system, arbitrary high error rates e can be tolerated. For 
mS e, the dynamics governed by the master equation 

Pt = m8£ s (p T ) + e£, no - iSC {p T ) 

is dominated by the first term and the steady state is 
arbitrarily close to the original steady state. In the spe- 
cific case, where the process Z s {pq-) — tr (pr)PT,s — P 
driving the target system into the steady state pq-, s is 
counteracted by depolarizing noise (tr(pj-)SL — pq-), the 
time evolution described by 

PT = mS(ti(pr)pr,s - p) +e(tr(p7-)II- p T ) 

leads to the steady state pfj- — m5 ^$l^ , which can be 
easily verified by solving the equation pj- = 0. This state 
is reached exponentially fast with a rate mS + e. The 
same result holds for local depolarizing noise acting on 
Alice's and Bob's system (see Sec. 4.1) if the steady state 
PT, S is a Werner state. A master equation of this type is 
solved exactly in the next section. 



4. SCHEME II: DISSIPATIVE ENTANGLEMENT 
DISTILLATION FOR WERNER STATES 



This is not generally the case, since the m source systems 
are coupled to each other through the target system. Due 
to this indirect coupling, the source systems may evolve 
differently in time and can reach different steady states, 
which can be disadvantageous for the evolution of the tar- 
get system. This is for example the case for the scheme 
described in Sec. 1.1 which does not include classical com- 
munication. 

It can be shown that ZsApt) — £s(pt)i if there 
is no state dependent back-action of T on the source 
systems. In this case, the evolution of the reduced 
density matrix of each source block is independent from 
the time evolution of the other blocks. This property can 
be guaranteed by re-initializing the source systems after 



In this section, we introduce a second dissipative dis- 
tillation scheme, which does not rely on entangling pro- 
cesses producing steady states, which are close to pure 
states, as scheme I presented in Sec. 1. We analyze here 
a very general model for Werner states [14], which can 
be solved exactly. Werner states are of the simple form 
Pw(f) = + (1 — /)(! — ^)/3, and are characterized 
in terms of their fidelity /, which is given by the over- 
lap with the maximally entangled state fi. Any quantum 
state can be transformed into a Werner state by twirling 
[14] without a loss of fidelity. Since a Werner-twirl is 
a LOCC map, a dissipative protocol can be constructed, 
which corresponds to the continuous application of a twirl 
operation on a given system and mapping of the resulting 
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state to a new pair acting as target system T by means 
of a continuous flip procedure (compare Sec. 4.2). 

This way, any dissipative process can be modified such 
that it can be described in terms of a Werner Lindblad 
operator Ef(p-j-) as used in Sees. 4 and 5, where / is the 
steady state fidelity of the underlying process. In this 
sense, the Werner model used here is very general and 
can be applied in many situations. 



4.1 Dissipative entangling model process for a single 
source pair 

A dissipative model process, which produces an ar- 
bitrary Werner state as steady state can be modelled 
by considering two processes, which generate the steady 
states SI and 1 respectively, where 1 =1/4 denotes the 
normalized identity. Let denote the four Bell-states, 
where | V'o) = ( |00) + |11)) /\/2, and cr, the Pauli matri- 
ces, where is the identity. A master equation which 
leads to the steady state £1 = |\&o) (^o| can be constructed 
using the four jump operators Qi = \ipo)(tpi\, which give 
rise to the Lindlbad term 

Q( P ) = J2£ Qi (p) = te(p)n-P- 

i 

Similarly, a master equation which leads to the steady 
statelis obtained using the jump operators Wij = (Xi®<Jj, 
which give rise to the Lindblad term 



W(p)=5> Wy (p)=tr(p)I- 



Hence, the Werner state pw(f) with fidelity / is the 
steady state of the time evolution governed by the master 
equation 

E f {p) = fQ(p) + ^l(3W(p) - Q( P )) = tx{p)p w - p. 

The Lindblad term Et(p) will be used in the following to 
model the basic entangling process acting on the source 
systems. 

Local depolarizing noise acting on Alice's (Bob's) 
side is included using the jump operators Si — Ia <£> &i 
(Si = Oi <£> Is), such that the corresponding Lindblad 
terms are given by 

N Mioe (p) = £ aL = PA <8Ib - P, 

i 

N Boh (p) = ^£f: b =lU®pB-/0, 

i 

where pa (Pb) is the reduced density matrix correspond- 
ing to Alice's (Bob's) system and I4 (H#) the normalized 
identity 1/2 on Alice's (Bob's) system. This process de- 
scribes the continuous replacement of the state on Alice's 



(Bob's) side by the completely mixed state. The total 
master equation 



i> = lE f {p) + -N{p) 



(15) 



where N(p) = N AUoe (p) + N Bob (p), describes the basic 
entangling process including local noise. This type 
of equation will be used frequently in the following 
sections, as it also describes also the evolution of the 
target systems once the corresponding source systems 
have reached the steady state. 

The steady state of the time evolution described 



by Eq. (151 is a Werner state 
7Pw(/) 



7 + e 

■y+s 



(16) 



The general time de- 



with reduced fidelity f s — 
pendent solution of the master equation ( 15 ) is of the 
form 

P{t) = Pogo(t) + Pi9i(t) + P292(t) + P39s(t), (17) 

where po is any initial state, pi = |(/c»o,a ®Ib+?a® Po,b), 
P2 = Pw(f) an d P3 =1 = Ia ®Ib- Po,a and Po.b are the 
reduced density matrices of the initial state po at Alices 
and Bobs side. The functions gi are given by 



9o(t) 
.91 W 

.93 (t) 



(< • 
1 



e 

2 I e 

1 

i 
1 



- e 

■i't 



-it 



(18) 



} (e- 1 ' 1 - e 7 "*) + ^(e 7 '* - + 1) 



where 7' = 7 + e and 7" = 7 + s/2. Note, that the terms 
which depend on the initial state of the system, i.e. po 
and pi, are suppressed exponentially fast. The system 



reaches the steady state given by Eq. ( 16 1 exponentially 
fast with a rate of at least 7 + f • 



In order to verify that Eqs. ( 17 ) and ( 18 1 are a solution 



of Eq. (15), Eq. (17 1 can be used as ansatz. The master 



equation gives rise to a set of differential equations for the 
functions gi with initial conditions go = 1 and gi — 0, 



9o 



■(7 + £)9o, 



(19) 



91 = -(7+ 2)51 + £ .9o, 

92 = -ef 2 +1(90 +9i +92 +.93), 



93 



-7.93 + 2-9i + £ .92- 



Below, the initial condition po = 1 will be considered 
frequently. In this case the solution simplifies to 



p(t) = p a + $-p s )e-^ t . 



(20) 
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4.2 Steady state entanglement distillation acting on 
n source systems 

We consider n systems which are subject to the basic 
entangling process jEf(p) + ^N(p) and are driven into 
the steady state p s as described in Sec. 4.1. These qubit 
pairs act as source systems for a LOCC distillation oper- 
ation To, which distills one potentially higher entangled 
state from these copies. The resulting quantum state is 
mapped to a target pair T and each source system is 
re- initialized in the state!. We do not specify To at this 
point - the solution derived in this section holds for any 
n to 1 distillation protocol. We start out by considering 
only deterministic protocols and generalize the results 
at the end of this section such that probabilistic schemes 
are also covered. Note that the complete re-initialization 
of the source systems represents the worst-case situation 
regarding the back-action of the target system onto 
the source pairs. This choice allows us solve the model 
exactly and to provide a lower bound for dissipative 
distillation schemes of this type. 

The continuous distillation procedure explained above is 
described by the master equation 



P : 



71 

E + l^frO). + 5 d(Td( P ) - p), (21) 



i=l 



where (X(p)) i stands for the dissipative process X(p) 
acting on the ith source system. 

In the following, we determine the time evolution 
and the steady state of the target system. The reduced 
master equation for T depends on the steady state 
of the reduced source system. Therefore, we start by 
solving the dynamics on the source system. Since the 
back-action of T on the source system does not depend 
on the quantum state of T, the time evolution of the 
source pairs can be considered independently from the 
target system. 

For clarity, the reduced states of source and target 
system are denoted by tr and pj- respectively in this 
section. The reduced master equation for the n source 
systems is given by 

n 

* = E + £ 2 N ^) l + Mtr(<x) I®" - a). (22) 

i— 1 

The solution of the homogeneous master equation which 
describes the entangling dynamics for n independent 
source systems 

n 

<7* (<J ,t) = E (l E f(< J *(< J Q, t )) +£ N {<7*(< J 0,t))) l , 
i=l 

is already known (see Sec. 4.1) if the initial state is a 
product state. <7* (ctq, t) denotes the solution of the homo- 
geneous master equation with initial condition er^oo, t — 



0) = (To- The solution of the inhomogeneous master equa- 
tion Eq. (22 1 is given by 

a(t) = a,(a ,t)e- SDt + So f dra,(f n , t - T )e- So( - t - r \ 

Jo 

= a*(a Q ,t)e- SDt + 5 D [ dra,(f n ,r) e - r ^ , 
Jo 

with arbitrary initial condition cr(0) = <7q. This solution 
can be easily verified by considering the time derivative 



&(t) = -6 D <7(t) + e- SDt <j*(a Q ,t) 
t 



So I dra* 
o 



, t - T)e 



t5t 



Using d t /J g(r)f(t - r) = f(0)g(t) + f* g(r)f(t - r) and 
a*(D®",0) =I g) ™, one obtains 

&(t) = -6 D a(t) + e- 6ot &*(a ,t) 

^0 



which yields Eq. (|22|. The steady state 



5n / dra* 



(23) 



is reached exponentially fast with a rate of at least So 
The homogeneous solution <j*(ff' n , r) is given by the ten 
sor product of the solution for a single source pair (20), 



a*(f n , t) = (p s + QL-p s ) e -^ t )® n , 



such that Eq. ( 23 1 can be further simplified 



dx(p s + (l-p s )x^)® n . 



(24) 



Next, we consider the dynamics of the target system T 
described by the time dependent master equation 

p T = S D (T D (cr(t)) - p T ) , 

which is solved by 

p r (t) = p T (0)e- 5ot + f drSoToia^e- 5 ^-^ 



with steady state Td(ct s ). The corresponding steady 
state fidelity can be inferred by integrating over the fideli- 
ties that are obtained if a standard distillation protocol 
is applied such that 

/out (f s ) = /out (/, e) = 1 1 dxfo (fs + ( 7 - Is )x ^ ) 8 " , 



where Eq. ( 24 ) was used. 



So far, it has been assumed, that the underlying 
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distillation protocol Tjj is deterministic, such that 
a distilled state is available whenever it is applied. 
However, many distillation protocols of interest are 
probabilistic, i.e., they only succeed some probability 
P(p). If a probabilistic distillation protocol is used, the 
corresponding map Tr> is defined in such a way, that a 
flip operation is only performed when the distillation 
was successful, which leads to a state dependent rate in 
the master equation 

p T = 6 D P(a{t))(T D (a(t))-p T ). 

Accordingly, the target system is driven into the same 
steady state as discussed above with a reduced rate. Once 
the time evolution of the source system has reached a 
steady state, the dynamics of the target system is deter- 
mined by the master equation 

p T = S D P(cr s ) (tr(p T )p' s - p T ) = 6 D pE fout (p T ), 

where p' s is the distilled steady state of the source system. 
Since p' s is a Werner state, the target system can act as 
one of n new source systems which drive a new target 
system into an even more entangled state. This way, the 
distillation protocol can be iterated in a nested form. 

5. CONTINUOUS QUANTUM REPEATERS 

The ability to distribute entangled states of high 
quality over long distances is of vital importance for 
quantum communication and quantum network related 
applications in general. As opposed to classical informa- 
tion, quantum information cannot be cloned. Therefore, 
classical repeater schemes are not applicable in this 
context and quantum repeater schemes which respect 
the coherence of quantum states are required [? ? ? 
]. In quantum repeater protocols, entanglement is first 
distributed over short distances L with high accuracy. 
Then neighboring pairs are connected by a teleportation 
procedure |S5 (entanglement swapping [? ? ]) such that 
entangled links which span a distance 2Lq are obtained. 
In the next step, two neighboring links of length 2Lq 
are connected by entanglement swapping, resulting in 
entangled pairs which span a distance AL$. This way, 
an entangled link of length L — Lo2 k can be established 
in k iteration steps (compare Fig. 4 in the main text). 
However, for non-maximally entangled states, entangle- 
ment swapping leads to a considerable degradation in 
the fidelity of the resulting quantum state. Since the 
distributed entanglement decreases dramatically every 
time the length of the entangling links is doubled, it 
can not be distributed over large distances this way. 
Therefore an entanglement distillation protocol has to 
be applied after every entanglement swapping procedure 
before proceeding to the next stage. 



t i r • 

m (fil) * 

Alice Bob Charlie 

FIG. S.6: (Color online) Continuous entanglement swapping. 

In the following, we describe a continuous dissipative 
quantum repeater scheme, which combines continuous 
swap and distillation processes in order to generate 
long-range entangled steady states, while entangling dis- 
sipative processes are only required over short distances. 
To this end, we introduce a continuous swap operation 
in Sec. 5.1 and explain in Sec. 5.2 how this method 
can be combined with the distillation scheme presented 
above (Sec. 4) such that a high-quality entangled link 
can be established over a large distance as steady state 
of a continuous dissipative evolution. We conclude this 
proof-of-principle study by giving a specific example. 



5.1 Continuous entanglement swapping 

The basic setup for entanglement swapping consists 
of three nodes aligned on a line, operated by Alice, 
Bob and Charlie, where Alice and Bob as well as Bob 
and Charlie share an entangled pair, while the distance 
between Alice and Charlie is too large for generating an 
entangled state of high quality (see Fig. 3b in the main 
text). By performing a teleportation procedure, which 
requires the measurement of the two qubits at Bob's 
node and classical communication to Alice and Charlie, 
as well as local operations on their sides, an entangled 
link can be established between Alice and Charlie [? ] . 

We consider a setting, where Alice and Bob as well as 
Bob and Charlie each hold a source pair which is subject 
to the basic dissipative entangling mechanism considered 
in Sec. 4, such both pairs are individually driven into the 
steady state p s . This dynamics is described by the Lind- 

blad term 7swE;=i ( E f(p))z = 7sw Ei=i( tr (pW - P)i- 
As illustrated in Fig. 3b in the main text, the source pairs 
are coupled to a pair of target qubits shared between 
Alice and Charlie through the term S sw (T sw (p) — p), 
where the completely positive map T sw corresponds to a 
flip operation which maps the state resulting from the 
entanglement swapping procedure to a target system 
and re-initializes the source systems in the state 1 <S> 1 
[S7] . Hence, the total master equation is given by 

2 

p = ^ E ( E M + l N 0>)) . + *«w (TsAp) - p) 

i=l 1 
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and the reduction to the source systems a yields 
2 

o = 7sw^ ( E f(p) + f^O )) - + *»w (tr(a)H®l- <r) . 

2=1 

The solution of this differential equation (compare 
Sec. 4.2) 



a{t) = a*(a ,t)e- 5 ™ t +6 l 



Jo 



where c^co^) is the homogeneous solution with initial 
condition a(t = 0) = oq, shows that the steady state 



ds(p, + 01- p a )a; ™w )« 



(25) 



j g reacnec i exponentially fast 



where p s = — 
with a rate of at least <5 SW . 

The time dependent master equation governing the dy- 
namics of the target system is given by 



P = S sw (T sw (a(t)) - p T ), 



(26) 



where T sw (a s ) is the steady state of this evolution. Ac- 
cording to Eq. (25), the steady state fidelity is given by 



f sw (f,e) = f sw (fs) = I dxF[f s + (]~ f s )x^ 



where f s is the fidelity of the state p s and F(f) — 
(1 - 2/ + 4/ 2 ) /3 is the output fidelity of the swap pro- 
tocol for two input states with fidelity /. A short calcu- 
lation shows that 

/sw(/i £) = /sw(/s) — 27 sw -t- -z — ^ r- 7(27) 

(27sw + o sw )(7 S w + o sw ) 

where / is the fidelity of the state pw and f s = 7s "^^ 4 ■ 
As discussed in Sec. 3, the scheme can be made robust 
against noise processes acting on the target system by 
using m copies of the source systems and coupling them 
all to the same target state. 



5.2 Creation of long-range, high-quality steady state 
entanglement 

The continuous swap operation introduced above 
(Sec. 5.1), the dissipativc distillation protocol explained 
in Sec. 4 and the method for stabilization against errors 
acting on target systems (Sec. 3) are the basic build- 
ing blocks for the dissipativc quantum repeater scheme 
illustrated in Fig. |S.7[ To begin with, the distance L 
over which an entangled link has to be established is di- 
vided into 2 k segments of length Lq, as in standard re- 
peater schemes. At each intermediate node, many qubits 



are supplied which are subject to local depolarizing noise 
acting at a rate e. We assume that each source pair con- 
stituting an elementary link of length Lq is individually 
driven into a steady state p s of high fidelity by means of 
an entangling dissipative process of the type discussed in 
Sec. 4.1., jEf T (p), with high initial steady state fidelity 
fi and a rate 7, which is large compared to the noise 
rate e. Note that this assumption can also be satisfied 
starting from dissipative processes leading to a steady 
state with low fidelity and low 7 if distillation and boost 
processes are applied as discussed above. In the follow- 
ing, we consider an iteration step of the repeater protocol 
which acts on 2 r entangled source systems, which each 
span a distance I with fidelity /; , and produces entangled 
links of the length 21 with fidelity f 2 i , such that f 2 i > fi ■ 
This is illustrated in Fig. |S.7| where the entangled source 
pairs of length I are shown in blue and the yellow target 
pairs of length 21 are depicted in yellow. Each iteration 
step consist out of the following subroutines, which are 
illustrated in Fig. [SJp: 



• Neighboring source pairs of length I (blue) are con- 
nected via a continuous swap operation. The re- 
sulting quantum states are written onto target pairs 
% Wt i (red) . In order to achieve a boost-effect on the 
targets, this protocol is run on m source systems in 
parallel. 

• A block of n such pairs %w,i,j , j = 1, . . • , n, acts 
as source system (green) for an distillation process, 
which maps the resulting quantum state to new 
target system 7em (yellow). 

• m of these blocks (green) are needed to achieve a 
high fidelity of the quantum state of the target sys- 
tems (yellow). 

This iteration step results in entangled links (yellow) 
which span twice the initial distance and feature a high 
fidelity as well as a high convergence rate once all source 
systems have reached the steady state. 

In the following, we consider 8d = S sw = ^- for 
simplicity (these parameters can be optimized for a 
given distillation protocol). The individual levels of 
the repeater scheme converge seriatim from bottom to 
top to a steady state. For example, once the source 
systems of length I (blue) are in a steady state, the 
reduced master equation for the target system of length 
21 (yellow) becomes time independent and this system 
reaches a steady state too. We assume now, that all 
source pairs of length I (blue) are driven by a time 
independent master equation of the type discussed in 
Sec. 4.2, p = jEf(p) once all underlying systems have 
reached the steady state. The reduced master equation 
for the target system T sw ,i (red) 



-PT„ 
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Lq Lq Lq Lq Lq Lq Lq Lq 



FIG. S.7: (Color online) Dissipative quantum repeater architecture, a) Concatenation of elementary steps, in which the distance 
over which entanglement is distributed is doubled, b) Illustration of a single iteration step including entanglement swapping 
and -distillation. 



includes local polarizing noise as introduced in Sec. 4.1 
and the back-action of the distillation scheme. Note, that 
the rate of the entangling process, Sm = 7 is again high, 
due to the boost on the target system. The entangling 
process acting on the target systems of the distillation 
procedure 7b , % (yellow), -E/ Tsw (pt ) is determined by the 
steady state fidelity /r 3W := f sw (fi,£)- 
A wide range distillation protocols for Werner states [S8] 
can be used in a continuous form as demonstrated in 
Sec. 4 (below a specific example is discussed). As ex- 
plained there, a distillation protocol corresponds to a 
completely positive map Td which is described by a Lin- 
blad term <5d (Td(/?) — p)- The distillation process is ap- 
plied continuously for each entangled link and the result- 
ing highly entangled qubit state is flipped to new target 
pairs 7b, z spanning the same length I. 
We consider here a distillation procedure which acts on 
n entangled source systems and distills one potentially 
higher entangled pair. Hence, for each of the 2 r_1 links, 
n copies % w> ij, i € {l,2 r ~ 1 }, j e {1,^} have to be sup- 
plied. This situation is sketched in Fig. |S.7b , where the 
target systems T sw .ij (red), driven by the source pairs 
(blue), are used as resource for creating a highly entan- 
gled steady state of the new target pair (yellow). One 
source block (shown in green) is sufficient for entangle- 
ment distillation, but several of them running in parallel 
are needed to boost the desired dynamics on the target 
system. This way, each target pair 7b, i is driven at a 
rate mfe = 7 and the total effective master equation for 
the target systems of the distillation protocol is given by 

PTo = mS D E fTo (p) + |iV(p). 

Hence, the resulting steady state fidelity is 

hu = (f D (f%J) = (/M/sw (//,£))) , 

where /d(/t sw ) is the entanglement distilled from source 
systems with steady state fidelity f%^- In order to iter- 
ate this process, we require fj- D > //, which can always 



be achieved using a strong entanglement distillation 
(large n) and high entangling rates 7 [S9] . 
The next iteration step begins with another continuous 
entangling swapping procedure. Here, the target systems 
of the distillation scheme act as source systems for the 
entanglement swapping operation. 

Since the total distance L has been divided into 2 fe 
segments ( L = L 2 k ), the protocol has to be iterated 
k times. As explained above, each iteration stage 
requires 2m 2 n qubit pairs such that in total (2m 2 n) 
source systems are needed. Hence the resources scale 
with (L/L ) log2(2m2,l) in the distance. We restrict 
the estimate of the required resources to the number 
of used qubit pairs, since the other resources scale 
polynomial in this quantity. As specific example we 
consider the distribution of an entangled state such that 
each repeater stage starts with and results in links with 
fidelity / = 0.96. We consider noise acing at a rate 
e = 0.05, distillation based on n — 16 source systems (the 
distillation protocol is described below) and stabilization 
of the target pairs by means of m — 50 copies of the 
underlying source blocks and 7 ~ 70. In this example, 
the required resources scale with (L/Lq) 16A - 

The entanglement distillation scheme used here is 
a four-to-one distillation protocol for Werner states |S8j 
which is applied two times in a nested fashion. Starting 
from four source states si, S2, S3, S4 with fidelity /i n , 
the following operations are performed. First, a bilateral 
CNOT Sl _>. S2 operation is applied to the first two pairs 
(where s\ is the control and S2 the target qubit ) and 
S2 is measured in the computational basis. Then, a 
Hadamard transformation is performed on both qubits 
of s%. Subsequently, a bilateral CNOT Sl _ ! . S3 operation 
is applied to the first and third pair (where si is the 
control and S3 the target qubit) and S3 is measured in 
the computational basis. The measurement obtained 
on Alice's and Bob's side are compared. If their mea- 
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surement results coincide, the resulting state Si is the 
desired higher entangled state. If not, the "safety-copy" 
S4 is used instead. In this event, distillation was not 
successful and the fidelity has not been increased, but in 
any case an entangled pair is available. The fidelity of 
the resulting state is given by 

, ,_ (l + g )(l + 7g 2 ) 

Jout \ Jin) 1 T P ' 

J-O-f; succ 

where g = (4/ <n - 1) /3 and P succ = (l + .g 2 + 2g 3 ) /4 is 
the success probability of this protocol. In the example 
above, this protocol is applied twice in a concatenated 
fashion. The second application of the scheme is run 
using the output states of the first one as input such 

that /£)(/,«) = font {font (fin))- 

We conclude this section with an estimate of the 
convergence speed of the presented repeater scheme. We 
start by considering the elementary pairs constituting 
the entangled links on the lowest level of the scheme. 
These systems reach the steady state up to a certain 
high accuracy after a time to. After this time, the 
dynamics of all systems on the next level is governed 
to a good approximation by a time independent master 
equation and converge with high accuracy to the steady 
state after another time period of length t has elapsed. 
Convergence of all k levels of the repeater scheme 
requires therefore a waiting time kt - Since the number 
of levels used in the scheme scale only logarithmical 
with the distance, we obtain a very moderate scaling 
of the convergence time with the distance. Note, that 
once the whole system is in a steady state, removal 
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